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Abstract 



We study a one-dimensional Frenkel Hamiltonian with off-diagonal dis- 
order, focusing our attention on the physical nature of the zero-energy peak 
of the density of states. The character of excitonic states (localized or delo- 
calized) is also examined in the vicinity of this peak by means of the inverse 
participation ratio. It is shown that the state being responsible for the peak 
is localized. A detailed comparison of the nearest-neighbor approach with the 
long-range dipole-dipole coupling is performed. 
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I. INTRODUCTION 



Since the pioneering works of Anderson,^ and Mott and Twose,il electronic and transport 
properties of randomly disordered systems have been the subject of long-lasting interest both 
from fundamental and applied viewpoints.^ - ! One-dimensional (ID) systems are frequently 
considered because they turn out to be simpler than those in three dimensions.il Originally, 
Mott and Twosei conjectured that all states are localized in ID systems, for any degree of 
disorder. Afterwards, a great deal of work has been devoted to examine Mott-Twose conjec- 
ture (see, for instance, the review |6|). However, it is well-known that electron derealization 
appears in ID random systems with short-range correlations.B! 

Two decades ago, Theodorou and Cohen established that the density of states (DOS) 
of a ID tight-binding Hamiltonian with nearest-neighbor (NN) interactions and random 
off-diagonal elements presents a singularity at the center of the band.! These authors used 
an analytical approach based on previous results obtained by Dyson0 for disordered linear 
chains of harmonic oscillators. In Ref. [|, it was also stated that the corresponding state is 
delocalized as the localization length was found to be infinite. Adding some amount of diag- 
onal disorder in the presence of off-diagonal randomness makes all states to be localized.!!!! 
Remarkably, the first calculations on ID tight-binding Hamiltonians with only diagonal dis- 
order did not reveal any singularity neither in the DOS nor in the localization lengthJillil 
Further, it was found both numerically^ and analyticallylllS a very weak anomaly (a peak 
but not a singularity) in both properties mentioned above. 

Recently, Fidder et al. have found by numerical diagonalization of the ID Frenkel Hamil- 
tonian with off-diagonal disorder that, notwithstanding the singularity of the DOS, the cor- 
responding state is localized if one includes the long-range (LR) interactions due to dipolar 
coupling between different sites.0 This finding seems to be in contradiction with the point 
of view raised in Ref. [5] suggesting that the state corresponding to the singularity of the 
DOS is delocalized. In this paper, we examine in detail the conclusions of Refs. |9] and [17. 
We address this issue by considering a ID Frenkel Hamiltonian with off-diagonal disorder 
with NN interactions and compare the results with those obtained when LR interactions are 
taken into account. The rest of the paper is organized as follows. In Sec. [Q], the ID Frenkel 
Hamiltonian with NN interactions is analyzed. We present arguments against those raised 
in Ref. [], namely that the zero-energy state is localized, even in the NN problem. This 
conclusion, based on analytical considerations, is then confirmed by direct diagonalization 
of the Hamiltonian. The detailed study of the ID Frenkel Hamiltonian with LR interactions 
is presented in Sec. fTJ. Section [TV] deals with the numerical simulations confirming the 
analytical results. Using numerical diagonalization of a Frenkel Hamiltonian with LR inter- 
actions, we calculate both the DOS and the inverse participation ratio, to be defined below, 
and study new features of these magnitudes with respect to the NN approach. Section [V| 
concludes with some comments regarding the results we have obtained. 



II. IS THE ZERO-ENERGY STATE DELOCALIZED? 

In this Section we briefly review the arguments of Ref. ^| leading to the conclusion that 
the state at the center of the band is delocalized. We present further arguments suggesting 
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the opposite point of view and, what is most important, numerical calculations confirm our 
statement. Let us consider a tight-binding Hamiltonian with only NN interactions 



H = Y, U n , n+1 (\n)(n + 1| + \n + l)(n\), (1) 

n 

where the NN interactions {U n;n+ i} are assumed to be 5-correlated and similarly distributed 
stochastic variables. The state vector \n) represents an excitation at site n. All site energies 
are set to zero since no diagonal disorder is included. The eigenvalue problem of the NN 
model reads 

Ea n , (2) 

where the set {a n } represents the real eigenvector corresponding to the eigenenergy E. For 
zero energy Eq. (||) gives the recurrence relation a n+ \ = — (C/ njn _i/C/ n>n+ i)a n _i. Using this 
relation one can find 

a 2 n+i= ( TJ "77 ••• [-JTr 1 - (3) 

\ ^2n,2n+l / \ ^ , 2n-2,2n-l / \ ^2,3 / 

The amplitudes at even positions vanish. The eigenvector (Q) represents the zero-energy 
state for a chain with odd number of sites. Defining the localization length at the center of 
the band L(E = 0) by the expression 



— lim — In 



L(E = 0) 2n 



a i 



(4) 



and applying the central-limit theorem, the authors of Ref. |9| obtained 1/L(E = 0) = 0. 
From this result they concluded that the state at center of the band was extended. 

The definition of the localization length (|]) is based on an unconditional assumption 
of the so-called exponential localization. Indeed, in such a case one would have typically 
«2n+i ~ exp[— (2n+l)/L]. Certainly, the definition (f|) cannot discern between a weaker than 
exponentially localized state (where the amplitude a2 n +i decreases or increases with n slower 
than an exponential) and an extended state (where d2 n +i ~ 1/ v^V with iV — > oo being the 
number of sites in the chain). In such a case, the mean extension of the eigenfunctions or 
the inverse participation ratio are the more adequate quantities for learning the character 
of the state. 

The results presented below show that the problem we are discussing just belongs to 
those which cannot be adequately analyzed from the assumptions leading to Eq. (|). First, 
let us write the NN interactions in the form U ntn +i = Uq(1 + £ n ,n+i), where £ n ,n+i is a 
Gaussian distributed stochastic variables with variance ^ < 1. Then, it is easy to calculate 
the probability distribution of ^2n+i = In |o-2n+i/ai|, 

£(Wi) = ^i-exp(-§±iY (5) 



From this, the authors of Ref. |TB| claimed that typically |o2n+i/ a i| ~ ex P(=t2£ 0v /rT). If 



so, one should conclude that the zero-energy state is localized rather than extended, in 
contradiction with the statement of Ref. |9|. 
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Further, using Eq. (|5|) we can calculate the probability distribution of X2n+i = |a2n+i/ai|, 

/(X2n+l) = r. ^ exp f- ln ,S" +1 ) • (6) 

This function has a sharp peak at Xmax — exp(— 2^n) and a very broad tail for large X2n+i 
such that (x2n+i) — I xf(x)dx — ex P(£o n )- Thus, it is rather difficult to make a definite 
conclusion from Eq. (^) concerning a typical dependence of |a2n+i/&i| on n. Nevertheless, 
the fact that /(0) = certainly indicates the zero probability to obtain an extended state. 
Below, we confirm this observation by numerical simulations. 



III. FRENKEL HAMILTONIAN 

We will be also interested in studying both the DOS and the degree of localization of 
states of a ID tight-binding Hamiltonian including all (LR) interactions, beyond the NN 
interactions. According to this, we then introduce the complete Hamiltonian 

N 

H = U mn \m)(n\, (7) 

m,n=l 

in which summation is performed now over all pair of sites. For definiteness, it is assumed 
hereafter that excitations described by the presented Hamiltonian correspond to Frenkel 
excitons. Furthermore, the U mn is assumed to be of dipole-dipole nature. We restrict 
ourselves to the case in which all transition dipole moments have the same magnitude and 
direction. Thus, one can take U mn = —U/\^ m — £ n | 3 , where —U (U > 0) is the dipole-dipole 
coupling of nearest-neighbors in the periodic lattice, i.e. at £ m — £ m +i = 1 (we chose here 
the negative sign of NN coupling as it takes place, for example, in J-aggregates0), and 
£ m = m + 5 m with 5 m being a stochastic variables assuming to be distributed around the 
regular sequence according to the Gaussian law with variance a 2 



A. The exciton spectrum and the DOS in the absence of disorder 

Before any discussions of the effects resulted from localization, it is useful to recall the 
peculiar features of the lD-exciton spectrum and of the DOS in the absence of disorder 
(5 m = 0). Then the Hamiltonian ([^D can be approximately diagonalized (with accuracy of 
the order of iV" 1 ) by introducing the excitonic basiJli 




The state vector \k) represents an exciton in the k-th state. Substituting @ into Eq. (^) 
one obtainsEl 
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N 

H = Y,E k \k)(k\, 

k=l 



(10a) 



E k = -2U± 1 cos (^) + O(N-i). (10b) 



Equation ( |10b| ) generalizes the corresponding expression of the NN approximation [n = 1 



term in ( |10b|) 1 to the case of including all (LR) interactions. We are especially interested 
in the behavior of the spectrum and of the DOS in the vicinity of extreme points, K = 
irk/ (N + 1) = and K = ir, as well as at the center of the band, k = (N + l)/2 (N taken 
to be odd). To do that, we exploit the following equation^ 

™ cos Kn f K\ 

£— = -H 2si,1 t) 

and the fact that one can extend the sum in Eq. (|10b|) up to infinity because n~ 3 decreases 
with n fast enough. Then, by integrating Eq. ( |TT| ) twice with respect to K, the sum in 
Eq. ( |10b|) can be cast into the formed 

E k = — 2Z7C(3) + -InlCjK 2 , K < 1, (12a) 



E k = 2^C(3) - U ln2(K - vr) 2 , K - n < 1, (12b) 

where ((3) = J2^Li n ~ z ~ 1.202. The corresponding formulae within the NN approximation 
are E k = -2U + UK 2 if K < 1 and E k = 2U - U{K - tt) 2 if K - tt < 1. Thus, one can 
conclude that LR interactions affect both the position of the bottom and the top of the band 
as well as the DOS of ID excitons. As we can see, the bottom and top of the band change, 
respectively, from -2U to -2E7C(3) » -2.40 4C/ and from 2U to (3/2)E/C(3) ~ 1.803C/. The 
DOS decreases approximately by the factor ^J\^l\E\ in the vicinity of the bottom of the band 
and, on the contrary, grows by the factor l/ln2 close the top. 

Finally, we would like to comment on the energy of the central exciton band state, with 
k = (N + l)/2. In the NN model, one finds that E( N+ iy 2 = 0, while with including all 
dipolar couplings, this energy is shifted to 

N 1 /7m\ 

E {N+1)/2 = ~2U £ — cos — ) « 0.225C/. (13) 

n=l n \ * / 

The DOS in the vicinity of the band center does not change noticeably as compared to the 
NN model. 
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B. Motion narrowing effect 



In the presence of disorder, the Hamiltonian of the system can be written as a sum of 
two parts: the unperturbed one ( |T0| ) and a term produced by the fluctuations of U mn 



N N 

H = Y,Ek\k)(k\+ J2 &kv\k)(k'\, 



(14a) 



k=l 



k,k'=l 



N 



^kk' 



N + 



7 E su v 



sm 



m,n=l 



irkn 
N+l 



sm 



nk'm 
N + l 



(14b) 



where SU mn = U mn — U where U means averaging over the probability distribution (||). 
Here A^r have diagonal and off-diagonal parts. The former is responsible only for the 
inhomogeneous broadening of excitonic levels, while the latter couples the excitonic modes 
and, therefore, causes the localization effects. 

The Afcfe/ undergo fluctuations because SU mn fluctuate. Assuming NN coupling and that 
SU mn fluctuations are small in some sense (see below), we can find the Akw distribution in 
an analytical form. This also helps us to comment on the results of numerical simulations 
that we discuss later in Sec. |V[ 

In order to achieve the task, we use the definition 



iV-l 



kk' 



5 I A kk > - 



N + 



j E w n , 



n+l 



n=l 



+ sin 




(15) 



Here, angle brackets indicate the average of the 5-function over the fluctuations of NN 
distances. They obey a Gaussian distribution law like (^) but replacing a 2 by a% N = 2a 2 . 
We omit the details of tedious but straightforward trigonometric calculations and only quote 
the final results. 

It can be easily shown that the sum in Eq. (|T5| ) is exactly equal to zero if k + k' = N + 1. 
Particularly, this means that A^k' = for k = k! = (N + l)/2 when iV is taken to be odd, 
i.e., the first order correction to the central energy is exactly equal to zero and does not 
fluctuate. Fluctuations of the other A k k' are distributed according to the Gaussian function 
with variances of the diagonal and off-diagonal elements distribution, aj(fc) and a 2 d (k,k'), 
given by 



(6aU) 2 
N + l 



2 + cos 



2-Kk 
N + l 



N + l 



(16a) 



(6aU) 2 
N+l 



1 + cos 



nk 



N + l 



cos 



' irk* ' 
N + l 



k + k' ^ N + l. 



(16b) 



From Eqs. ( |T6"D one can conclude that, in the case of off-diagonal disorder, the motion 
narrowing effect is also present as it takes place for diagonal disorder,!^ i.e., both magnitudes 
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in (|16D scale as (N + I)' 1 - We should point out that, in contrast to diagonal disorder, here 
the magnitudes a d and a nd are functions of the state numbers. Note also that <Td(k) goes 
through its minimum value exactly at the center of the exciton band, i.e., at k = (N + l)/2 
and k = N/2 at N taken odd and even, respectively. In fact, we can also assert this with 
respect to the value of a n d(k, k') since k and k! cannot differ greatly provided the condition 
a nd (k, k') < U. 

To conclude this Section let us comment on the validity of the motion narrowing effects. 
Obviously, this is valid only when a nd < \Ef. — E^ + i\. In this case, the excitonic states 
are not mixed by the perturbation and remain to be extended over the whole chain. They 
are essentially mixed for the opposite sign of inequality, reducing their localization lengths. 
Then, the number of sites within the region of localization (N*) drives the motion narrowing 
effect rather than the whole number in the chain N. In Refs. |TJ] and |22] a self-consistent rule 
for estimation of N* is carried out. 



IV. NUMERICAL SIMULATIONS AND DISCUSSIONS 

Further, we will mainly focus our attention on the normalized density of states p(E) and 
on the degree of localization (inverse participation ratio, IPR) for the states at energy E. 
They are defined respectively as followsS 




(17a) 



where the angular brackets indicate an average over an ensemble of disordered linear chains 
and the a^ n is the eigenvector of the Hamiltonian (0) corresponding to the eigenvalue Ej, 
with k = 1,2,3, ...,iV, 

N 

U nm a km = E k cikn- (18) 

m=l 

The IPR behaves like 1/N for delocalized states spreading uniformly over the entire system 
on increasing N. In particular, the IPR can be exactly computed for the eigenstates of the 
periodic lattices given in (|9]). In doing so we obtain the expected behavior for N —>■ oo. 
On the contrary, localized states exhibit much higher values. In the extreme case, when the 
exciton is localized at a single site, the IPR becomes unity. Therefore, the scaling analysis of 
the IPR as a function of the system size provides valuable information about the nature of 
the excitonic states. We should mention that a complete multifractal analysis, accomplished 
by studying the scaling of the other moments of the probability distribution, is beyond the 
scope of this work. 
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FIGURES 



1 1 
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c) 





-4 -2 2 4 
E 

FIG. 1. Density of states in the frame of NN coupling when the lattice size is N = 2500 and 
the degree of disorder is a) &NN = 0.02, b) 0.04, and c) 0.16. 

We have solved numerically the eigenvalue problem (|I8|) for different values of disorder, 
namely the mean fluctuation of the NN distance, <jnn = V%cr, to study the features of 
both the DOS and the IPR discussed above. In our numerical treatment <jnn ranges from 
(periodic lattices) up to 0.32 whereas the maximum system size we have considered is 
N = 2500. Results comprise averages over 50 realizations of the disorder for each given pair 
of parameters N and (Jnn- 



A. Nearest-neighbor approximation 

Let us comment the results we have obtained for the NN approximation. Figure [I] shows 
the DOS for the largest lattice size we have considered (N = 2500) and different values 
of the disorder (cnn = 0.02,0.04,0.16 from top to bottom). We observe that the DOS 
is symmetric about the center of the band. The singularities at the edge of the exciton 
band are smeared out on increasing the degree of disorder. Interestingly, a sharp peak in 
the DOS at the center of the exciton band appears when the degree of disorder exceeds 
some threshold value (<tnn ~ 0.03 for our model parameters). We will discuss further this 
point later. We have also observed that the percentage of states in the DOS peak increases 
with the degree of disorder. In addition, the amplitude of the peak rises noticeably with 
increasing the number of sites in the lattice, as seen in Fig. ||] for ct^n = 0.08. 
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FIG. 2. Density of states in the frame of NN coupling when the degree of disorder is &nn 
and the lattice lattice size is a) N = 1000, b) 1500, and c) 2500. 
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The IPR presents an overall increase when the degree of disorder increases, meaning that 
the larger the degree of disorder, the smaller the exciton localization length. This is clearly 
observed in Fig. |3|, where we show the IPR as a function of the exciton energy for the same 
parameters of Fig. [1]. However, the increase of the IPR strongly depends on the energy, being 
more pronounced close to the center of the band. Simultaneously with the occurrence of the 
peak of the DOS, a hardly visible downfall arises in the IPR at zero energy. This downfall 
is better revealed for larger lattices, as it can be seen in Fig. [| for the same parameters of 
Fig. §. 




FIG. 3. Inverse participation ratio for the same cases shown in Fig. |]. Notice the overall 
increase on increasing the degree of disorder. 
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FIG. 4. Inverse participation ratio for the same cases shown in Fig. [|. The inset shows an 
enlarged view of the center of the band. 



As mentioned above, the scaling of the IPR with the lattice size may be useful to discern 
the nature of the eigenstates. The IPR at the center of the band for different values of 
the degree of disorder is plotted in Fig. [5] as a function of the lattice size. The IPR for 
periodic lattices scales very accurately as 1/N, hence indicating that their eigenstates spread 
uniformly over the whole lattice. As soon as some amount of disorder is introduced in the 
system, the IPR follows a power law for small N but tends to a constant value for large N, 
as plotted in Fig. |5|. The critical size for which deviation from power fit occurs decreases 
on increasing the degree the disorder. The constant value of the IPR for large N increases 
with the degree of disorder, indicating that the eigenstates at the center of the band actually 
become more localized. 
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FIG. 5. Scaling of the inverse participation ratio with system size for the eigenstates at the 
center of the band. Labels indicate the degree of disorder. 

Summarizing these observations for the NN approximation, we are led to two main 
conclusions. First, the zero-energy peak of the DOS really exists. Moreover, as its width 
shows no dependence on the degree of disorder (at least, when the latter ranges over the 
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interval used in our simulations), we are inclined to identify this peak with a 5-like singularity 
rather than to the famous Dyson singularity ~ 1 / 1 In 3 \E\. This singularity was found 
firstly for a special distribution of the NN hopping integral in the form of a generalized 
Poisson function.0 Second, the corresponding eigenstates show no tendency to derealization 
with rising the lattice size contrary to the opposite statement done in Ref. ^. Moreover, they 
are not more delocalized than those of energies close to zero. 



B. All interactions 



Effects of inclusion of all dipolar interactions in Eq. (|18D on the DOS and IPR has been 
already discussed in Ref. [l?]. Nevertheless, it has been done only for a fixed values of the 
chain length (N = 250) and the degree of disorder (a = 0.08). Below we present our DOS 
and IPR numerical data obtained by varying both N and a. 

Figures ^] and [?] show the results of numerical calculations of both the DOS and the IPR 
for different values of the degree of disorder and N = 2500. Here, one can observe the usual 
changes of both magnitudes as compared to those in the NN approximation: asymmetry 
and shift of the excitonic band edges, both in a good agreement with the analytical results 
presented in Sec. |111 A . 

In addition, some new features appear, namely the peak in the DOS has a finite width 
and is shifted from zero-energy to a somewhat higher value -E pea k ~ 0.21U for low degree 
of disorder, in full correspondence with the results of numerical simulation done in Ref. [T7| . 
Higher values of the degree of disorder lead to a smaller energy shift. Further, a peak in 
the IPR appears at the same energy where the DOS peak, with a finite width as well. The 
last observation means that the states forming the DOS peak become more localized as 
compared to those with close energies, in contrast to the case of the NN interaction. This 
is also confirmed by the scaling of the IPR for -E pea k with the system size (not shown here): 
In all cases we observe higher values of the IPR in comparison with those shown in Fig. ||. 
Besides that, the trend is similar, that is, the IPR scales as 1/N only for perfect lattices, 
whereas it tends to a constant value for non-zero degree of disorder. 
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FIG. 6. Density of states in the frame of LR coupling when the lattice size is N = 2500 and 
the degree of disorder is a) onn = 0.02, b) 0.04, and c) 0.16. 
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FIG. 7. Inverse participation ratio for the same cases shown in Fig. ||. 

C. Discussion 

Now let us discuss the origin of the features of the DOS and of the IPR found in numerical 
simulations. 

1. NN interaction 



Obviously, the zero-energy peak in the DOS might appear when the states at the center 
of the exciton band become localized, i.e., their localization lengths are reduced to values 
less than the lattice size. This occurs when the reduced degree of disorder due to motion 
narrowing, (7^(0,0) = Qo~nnU / (N + l) 1 / 2 , exceeds the energy spacing at the center of the 
band, AE = 2ttU / (N + 1). Equalizing these two magnitudes one obtains an estimation for a 
threshold of mean fluctuations of the NN distances to observe the peak, a^} N ~ U/(N+ l) 1 / 2 . 
Thus o~j!} N ~ 0.02 for iV = 2500, which is in a good agreement with the numerical data of 
Fig.0. 

With regard to the fact why this peak appears, we can suggest two explanations which 
seem to be suitable for the model under consideration. First, as the distribution of disorder 
we used has long tails then, owing to possible large fluctuations of the NN distances, strongly 
interacting dimers can be created whose level splittings noticeably exceed the typical mag- 
nitude of the intersite interaction U. Consequently, the whole chain is broken into several 
independent segments in the sense that two adjacent dimers produce a potential well for the 
exciton, localizing it into the segment bounded by them. As the zero eigenenergy is always 
present in a segment with odd number of sites, one can expect a peak in the DOS at this 
energy (similar explanation of the zero-energy peak of the DOS was suggested in Ref. |18| ). 
The peak amplitude increases with disorder simply because of the rising of the number of 
segments as the degree of disorder grows. Appearance of such strongly interacting dimers is 
clearly seen from that the IPR approaches 0.5 at the DOS tails (see Figs. |3] and |]). 

Recently, it was demonstrated that the Dyson singularity of the DOS appeared even for 
a box-like distribution of disorder.il Then, the explanation above fails due to the absence 
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of large fluctuations of the NN-randomness at a low magnitude of the degree of disorder. 
In such a case, another cause for the occurrence of the zero-energy peak in the DOS can 
be proposed. As we have already noted in subsection [Tl [JJ, the first order correction to the 
central energy is equal to zero for chains with odd number of sites and has a minimum of 
fluctuation in the case of even number of sites. The zero-energy peak indicates that the 
central band eigenenergies are more stable to perturbations than the remaining ones. This 
certainly will result in a peak of the DOS after averaging over realizations of disorder. It is 
remarkable that simulations done for a special type of disorder — which has no effect on a 
certain excitonic level in the sense that the first order correction to the energy vanishes — 
show an analogous peak in the DOS at this energy.fEl Thus, this empirical rule can serve for 
inspecting the appearance of peaks in the DOS for the tight-binding Hamiltonian. As the 
last treatment does not use any specific peculiarities of the NN-randomness distribution, it 
seems to be suitable for any other distribution. We suppose that, for the model considered in 
this paper, both mechanisms discussed above contribute to the formation of the zero-energy 
peak in the DOS. 

Concluding this subsection, note that the degree of localization of the central states 
obtained from the numerical simulation is in a good agreement with the theoretical estimates 
based on a self-consistent rule proposed in Refs. [TP] and |22] (see paragraph 1 of the present 
subsection) . 



2. All interactions 

As it was stated in Ref. [T7|, the energy shift of the DOS peak, E veak ~ 0.21C7, agrees 
very well with the energy of the central band state in the absence of disorder, Eq. (|T3|). 
We are also inclined to relate this peculiarity to a state of analogous origin, i.e., similar to 
sin(7m/2). This can be demonstrated at least in the perturbative limit. Moreover, exploiting 
this analogy further, we should assume that the character (having no amplitude on the half of 
sites) of the mentioned eigenstate has not to be changed dramatically (at least, in average), 
when passing from the NN model to the exact one, as it is the case for the problem without 



disorder.! 
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The 5-singularity of the DOS becomes broader with including all dipolar couplings. At 
least two effects can contribute to this broadening. As it was supposed above, the DOS 
peak results from the isolated segments of odd number of sites, which, in turn, originate 
from large fluctuations of the NN distances. At moderate magnitudes of disorder we are 
mainly dealing with, it is unlikely the simultaneous strong reduction of the distance between 
a nearest-neighbor pair and the distances with other neighbors. Thus, for the very beginning, 
one can consider adjacent segments as independent of each other. Then, the eigenenergy 
of the local (belonging to a certain segment) central band state will fluctuate, owing to 
fluctuations of the segment lengths [see Eq. (pT3|)1, and thus will produce inhomogeneous 
broadening of the DOS peak. The second probable origin of this effect is the coupling of 
different isolated segments due to the interaction with far neighbors. 



In Ref. [L7|, the appearance of the IPR peak was explained by an exceptional property 
of this characteristic with regard to the central band state, k = (N + l)/2, characterized 
by the wave function [2/ (AT + l)] 1 / 2 sin(7m/2). Even in the absence of disorder, the IPR 
defined by Eq. ( p. 7 bp shows stronger localization of this state [£ = 2/(N + 1)] as compared 



13 



to localization of the remaining states [C = 3/2 (AT + 1)].0 The authors of Ref. [TJ] asserted 
that the IPR peak in the presence of disorder reflected a remnant of this special state in 
those forming the peak. At the moment, we do not see any other explanation of the origin 
of this anomaly. If it is so, the similar feature might be manifested in the NN problem, too. 
Nevertheless, as follows from our simulations done for the NN problem, the IPR displays a 
downfall rather than a peak. One of the reason for such difference may be the fact that the 
zero-energy state is not exponentially localized in the NN problem (see Sec. [TT]). It results in 
a larger extension of this state as compared to the others. In principle, such a large extension 
can compensate the IPR anomaly coming from the special character of the zero-energy state 
(having no amplitude at all on the half of sites) giving rise to the same value of the IPR at 
E = and at surrounding energies. 

V. SUMMARY 

The numerical study of the problem of the zero-energy peak of the DOS for a one- 
dimensional Frenkel chain with only off-diagonal randomness shows that the peak is really 
presented. In the NN approximation, it is located at the center of the excitonic band and 
tends to convert to a (^-singularity as the size of the chain increases. The states belonging 
the peak are localized and do not display any tendency to derealization with the chain 
size. Moreover, the degree of localization (IPR) does not differ very much from that of the 
surrounding states. The inclusion of couplings due to far neighbors shifts the peak to a 
slightly higher energy 0.21U), while the IPR, in contrast to the NN problem, shows a 
peak at the same energy. 
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